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Abstract. The production of mesons in proton-nucleus collisions from 1.0 to 2.5 GeV is analyzed 

' with respect to one-step nucleon-nucleon (7VA'' NYK^) and two-step Zi-nucleon {AN — > K'^YN) or 

pion-nucleon {tvN K^Y) production channels on the basis of a coupled-channel transport approach 

, (CBUU) including the kaon final state interactions. The infiuence of momentum-dependent potentials 

for the nucleon, hyperon and kaon in the final state are studied as well as the importance of K'^ elastic 
rescattering in the target nucleus. The transport calculations are compared to the experimental K'^ spectra 

! taken at LBL Berkeley, SATURNE, CELSIUS, GSI and COSY- Jiilich. It is found that the momentum- 

dependent baryon potentials effect the excitation function of the cross section; at low bombarding 
energies of ~ 1.0 GeV the attractive baryon potentials in the final state lead to a relative enhancement 
of the kaon yield whereas the net repulsive potential at bombarding energies ~ 2 GeV causes a decrease 

' of the cross section. Furthermore it is pointed out, that especially the spectra at low momenta 

(or kinetic energy Tk) allow to determine the in- medium potential almost model independently due 

, to a relative shift of the K'^ spectra in kinetic energy that arises from the acceleration of the kaons when 

propagating out of the nuclear medium to free space, i.e. converting potential energy to kinetic energy of 
the free kaon. 

PACS. 13.60.Le Meson production - 13.75.Jz Kaon-baryon interactions - 14.40.Aq Pi, K, and eta mesons 
- 24.40.-h Nucleon-induced reactions 



1 Introduction 

The production of mesons heavier than pions in p + A re- 
actions at bombarding energies far below and close to the 
free nucleon-nucleon threshold is of specific interest 0]- 
]l8t as one hopes to learn either about cooperative nuclear 
phenomena and/or about high momentum components 
of the nuclear many-body wave function that arise from 
nucleon-nucleon correlations. Especially mesons have 
been considered as promising hadronic probes |19,2^ due 
to the rather moderate final state interaction, which is a 
consequence of strangeness conservation and the fact that 
there are no baryon resonances with anti-strange quarks 
in nuclei. Anti-hyperons, furthermore, have a much larger 
production threshold and annihilate very fast in nuclei. 
On the other hand, the kaon properties might change in 
the nuclear medium such that conclusions on cooper- 
ative nuclear phenomena require a precise understanding 
of the (anti-) kaon potentials at finite nuclear density. 

Experiments on production from nucleus-nucleus 
collisions at SIS energies of 1~ 2 A- GeV have shown that 
in-medium properties of the kaons are seen in the collec- 
tive flow pattern of K'^ mesons both, in-plane and out- 



of-plane, as well as in the abundancy of antikaons p^,p3[. 
Thus in-medium modifications of the mesons have become 
a topic of substantial interest in the last decade triggered 
in part by the early suggestion of Brown and Rho |24| , that 
the modifications of hadron masses should scale with the 
scalar quark condensate < qq > at finite baryon density. 

As demonstrated in the pioneering work of Kaplan and 
Nelson kaons and antikaons couple attractively to the 
scalar nucleon density with a strength proportional to the 
KN — S constant, 

Skn ^]^{ml+ml) <N\uu + ss\N>, (1) 

which is not well known at present and may vary from 270 
to 450 MeV. In (|]) to^ and m° denote the bare masses for 
the light u- and strange s-quark while the expression in 
the brackets is the expectation value of the scalar light 
and strange quark condensate for the nucleon |Q. Fur- 
thermore, a vector coupling to the quark 4-current - for 
vanishing spatial components - leads to a repulsive poten- 
tial term for the kaons; on the other hand this (Weinberg- 
Tomozawa) term is attractive for the antikaons. 
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In chiral effective tlieories tfie dispersion relation for 
kaons and antikaons in the nuclear medium - for low mo- 
menta - can be written as 125(1 



WK(pAr,p) = +C, ^2 
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In equations (2), (3) ttik denotes the bare kaon mass, 
« 93 MeV is the pion decay constant, while Ps and p^ 
stand for the scalar and vector nucleoli densities, respec- 
tively. As shown in Ref. for Ekn= 450 MeV one ends 
up with an effective kaon potential which is close to zero 
below ordinary nuclear matter density po and becomes 
more repulsive above po- On the other hand, using IJkn= 
270 MeV a repulsive kaon potential of w 25 MeV at nor- 
mal nuclear matter density is obtained. Note, that when 
extrapolating to 3po and above, the antikaon mass be- 
comes very light. Thus antikaon condensates might occur 
at high baryon density which, furthermore, are of great 
interest in the astrophysical context p7|, p8||29[ |. 

However, the actual kaon and antikaon self energies (or 
potentials) are quite a matter of debate ~ due to higher 
order terms in the chiral expansion - especially for the an- 
tikaon |^,^|3^ and the momentum-dependence of their 
self energies is widely unknown (except for a dispersion 
analysis in Ref. |Q) since most Lagrangian models re- 
strict to s-wave interactions or only include additional p- 
waves. It is thus mandatory to perform experimental stud- 
ies of the (anti-) kaon properties under well controlled con- 
ditions, e.g. in proton-nucleus reactions, where one probes 
the (anti-) kaon self energies at normal nuclear matter den- 
sity po !v 0.16 fm^"^ and below. Furthermore, by gating on 
kaon momenta in the laboratory, one might be able to ob- 
tain information on the momentum dependence of the self 
energies, too. 

production m p + A collisions at subthreshold 
energies has been observed experimentally more than a 
decade ago by Koptev et al. |20) at bombarding energies 
from 0.8-1.0 GeV. However, only total K'^ yields could 
be extracted at that time. Nevertheless, the target-mass 
dependence of the yield (~ A) suggested the domi- 
nance of two-step reactions with an intermedia te p ion or 
A. Detailed folding-model calculations in Refs. |1^,M es- 
sentially came to the same conclusion. First differential 
spectra from p + NaF and p + Pb reactions from 
the LBL Berkeley had been performed at Tiab = 2.1 GeV 
g, i.e. far above the NN threshold of 1.58 GeV in free 
space. Only in more recent years differential spectra 
have been measured down to 1.2 GeV for ^^C targets at 
40 ° H (at SATURNE) or 90 ° in the laboratory || 
(at CELSIUS). Unfortunately, the different experiments 



have no overlap in acceptance and the interpretation of 
the data, if compatible at all, remains vague ||l^. First 
data on the full momentum distribution at forward angles 
have been presented very recently by the ANKE Collabo- 
ration at COSY-Jiilich |37| for mesons from p +^'^ C 
reactions at 1.0 GeV [|38| . These data show a kinemati- 
cal focussing of the spectra at finite momentum of ~ 350 
MeV/c, which appears incompatible with the cross section 
from Ref. . Thus a systematic comparison of all data is 
urgently needed within an adequate theoretical approach 
that allows to compare the kinematically restricted data 
on the same footing. 

In this study we use the coupled-channel (CBUU) trans- 
port model that has first been developed in Ref. [g9| for 
the description of nucleus- nucleus collisions and later on 
employed for the simulation of pion- and proton-nucleus 
reactions [^,40,4l|, too. For applications to produc- 
tion in nucleus- nucleus collisions at SIS energies we re- 
fer the reader to Ref. In this model the effects of 
momentum-dependent self energies for all hadrons can be 
studied explicitly as well as their production and propa- 
gation in the nuclear medium. 

The paper is organized as follows: We briefiy recapit- 
ulate the ingredients of the CBUU model in Section 2, 
present the extensions performed in this study and in- 
vestigate in particular the effects from rescattering. 
In section 3 we compare the transport calculations to 
the experimental spectra available from different labora- 
tories and explore the sensitivity of the K'^ spectra to the 
momentum-dependent potentials employed. A summary 
and discussion of open problems concludes this paper in 
Section 4. 



2 Ingredients of the transport model 

In this work we perform our analysis along the line of the 
CBUUQ approach |39| which is based on a coupled set 
of transport equations for the phase-space distributions 
fh{.x,p) of hadron h, i.e. |||,|4) 



+ {n^d^u'^, + Mid^ui) d^Ux,p) 
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~ fh{x,p)fh2{x,P2)fh3{x,P3)ffu{x,P4)} ■■■ ■ (5) 

In Eq. (^) U^(x,p) and Ul^{x,p) denote the real part of the 
scalar and vector hadron selfenergies, respectively, while 
[G''G]i2^34.J^{n + 772 - 7^3 - 7/^4 . . .) is the 'transition 
rate' for the process 1 + 2^3 + 4 + ... which is taken to 
be on-shell in the semiclassical limit adopted. The hadron 
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quasi-particle properties in 
shell constraint 



are defined via the mass- 
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with effective masses and momenta (in local Thomas- Fermi 
approximation) given by |Q 



n^{x,p) = p- 

while the phase-space factors 



A'h + U^{x,p) 

'^-ui:{x,p) 



fh{x,p) = 1 ± fh{x,p) 
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are responsible for fermion Pauli-blocking or Bose enhance- 
ment, respectively, depending on the type of hadron in the 
final/initial channel. The dots in Eq. (5) stand for further 
contributions to the collision term with more than two 
hadrons in the final/initial channels (cf. Ref. |45)). The 
transport approach (^) is fully specified by U^{x,p) and 
Uj^{x,p) (/i — 0,1,2,3), which determine the mean-field 
propagation of the hadrons, and by the transition rates 
G''G6'^{...) in the collision term (5), that describes the 
scattering and hadron production/absorption rates. 

The scalar and vector mean fields and for nu- 
cleons are modeled as in Ref. , however, slightly mod- 
ified in line with Ref. . In Fig. |l] the real part of the 
Schroedinger equivalent potential (SEP) for nuclcons. 



Use pip) - ir°(p) - x/p- 
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is shown at density po (full line) as a function of the mo- 
mentum p relative to the nuclear matter at rest. Whereas 
we see a net attraction for momenta p < 0.6 GeV/c, the 
nucleon potential becomes repulsive for higher momenta 
and reaches a maximum repulsion at p ~ 1.5 GeV/c. We 
mention that at density po the Schroedinger equivalent 
potential Usep compares well with the potential from 
the data analysis of Hama et al. |^ as well as Dirac- 
Brueckner computations from up to a kinetic energy 
Ek^n of 1 GeV @. 

Apart from the nuclear potentials each charged hadron 
additionally moves in the background of the Coulomb po- 
tential that is generated by the charged hadrons them- 
selves. In case of proton- nucleus reactions - with the nu- 
cleus at rest ~ this essentially amounts to a Coulomb ac- 
celeration in the final state for positively charged hadrons. 
Note, that for heavy nuclei like Pb or Au the Coulomb po- 
tential in the nuclear interior is about -I- 20 MeV, i.e. of 
the same order of magnitude as the 'expected' repulsive 
nuclear potential. 

The hyperon mean fields, furthermore, are assumed 
to be 2/3 of the nucleon potentials. In the present ap- 
proach, apart from nucleons, Zi,iV(1440), N{1535), A,S 
with their isospin degrees of freedom, we propagate ex- 
plicitly pions, K~^,K~, and ry's and assume that the pi- 
ons as Goldstone bosons do not change their properties in 
the medium; we also discard self energies for the 77-mesons 
which have a minor effect on the kaon dynamics. The kaon 
and antikaon potentials, however, have to be specified ex- 
plicitly. 
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Fig. 1. The nucleon potential (solid line) at density po as a 
function of momentum relative to the nuclear matter at rest 
as used in the CBUU transport approach. The hatched area 
denotes the K'^ potential at po (within the uncertainties) from 
the dispersion analysis in Ref. |33]. 



2.1 K^,K self energies 

Apart from (||) , ( 3|) there is a couple of models for the kaon 
and antikaon se fenergies |2l|,^,|49|, which differ in the 
actual magnitude, however, agree on the relative signs for 
kaons and antikaons. Thus in line with the kaon- nucleon 
scattering amplitude the K'^ potential should be slightly 
repulsive at finite baryon density whereas the antikaon 
should see an attractive potential in the nuclear medium. 
Without going into a detailed discussion of the various 
models (cf. Ref. and figures therein) we adopt the 
more practical point of view, that the actual and K~ 
self energies are unknown and as a guide for our analysis 
use a linear extrapolation of the form, 
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with a{p) denoting a momentum-dependent coefficient. 
Following the dispersion analysis of Sibirtsev et al. [|3| 
the coefficients a{p) can be modelled in line with the 
potentials from Fig. 9 of |Q; the resulting kaon poten- 
tial Ux+ at density po is shown in Fig. |^ in terms of the 
hatched area and remains repulsive for all momenta con- 
sidered. In the momentum range of interest in this study, 
i.e. from 0.1 ~ 1.0 GeV/c, the potential at density 
Po may be approximately represented by a constant of 
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Uk+ ~ 20 MeV taking into account the relative uncer- 
tainty of ± 10 MeV from the analysis in Ref. |33|. Since 
the antikaon dynamics has been investigated in Ref. p3| 
for p + A reactions, we skip a further description of the 
actual implementation of the K~ potential. 



2.2 Perturbative treatment of strangeness production 

The calculation of 'subthreshold' particle production has 
to be treated perturbatively in the energy regime of in- 
terest due to the small cross sections involved. Since we 
work within the parallel ensemble algorithm, each parallel 
run of the transport calculation can be considered approx- 
imately as an individual reaction event, where binary re- 
actions in the entrance channel at given invariant energy 
y/s lead to final states with 2 (e.g. K~^Y in ttB channels), 
3 (e.g. for K~^YN channels in BB collisions) or 4 parti- 
cles (e.g. KKNN in BB collisions) with a relative weight 
Pi for each event i which is defined by the ratio of the 
production cross section to the total hadron-hadron cross 
section^ We thus dynamically gate on all events where a 
K^Y or K^K~ pair is produced initially. Each strange 
hadron production event i is represented by ~ 10^ test- 
particles for the final strange hadron j with individual 
weight Wj such that the sum of the weights over j re- 
produces the individual production probability Pi and the 
distribution in momenta (multiplied with the A^iV or ttN 
cross section) describes the differential production cross 
section. Then the 'perturbative' hadrons are propagated 
according to the Hamilton equations of motion including 
the potentials. Elastic and inelastic reactions with pions, 
?7's or nonstrange baryons are computed in the standard 
way |§. 

The final differential cross sections are obtained by 
multiplying each testparticle weight Wlj by the total in- 
elastic pA cross section and gating on the experimen- 
tal acceptance of the different detectors. In this way one 
achieves a realistic simulation of the strangeness produc- 
tion, propagation and reabsorption during the proton-nucleus 
collision with sufficient statistics to allow for selective cuts -'^'S 
also at the low bombarding energy of 1.0 GeV, where the 
total X"*" cross section is in the order of 1 fxh (for Au) or 
below (for 12 C). 



correspond to the experimental data for the different ttN 
channels in vacuum (or 'free' space). Note, that in addi- 
tion to our early studies in |]13j[i1] the channels with a 
S hyperon are taken into account. All cross sections are 
reparametrized as a function of the invariant energy above 
threshold ^/s — ^/sq , where y/s^ denotes the threshold 
for the individual channel given by the sum of the hadron 
masses in the final state of the reaction. 

Whereas the production cross section of kaons from 
pN collisions close to threshold has been essentially un- 
known about a decade ago, the ambiguity in this cross 
section has been resolved experimentally by now and 
more adequate parametrizations of the cross sectio ns can 
be employed. The experimental data from Refs. 
on the pp + X reactions are displayed in Fig. 

comparison to the current approximation (solid line) and 
the estimate from Zwermann (dotted line) used earlier 
[|7[|l|. Thus the problem of 'free' cross sections for pN 
collisions is sufficiently under control. For AN collisions, 
however, no experimental data are available. We use the 
same cross section as a function of the invariant energy ^/s 
as in the pN case keeping in mind this basic uncertainty. 
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2. Comparison of our parametrization for the inclusive 
production cross section in pp reactions (solid line) with 
the data from Refs. [^,^. The dotted line corresponds to the 
parametrization from Zwermann used earlier in Refs. 



2.3 Elementary cross sections 

For the present study the production of pions by pN col- 
lisions in p + A reactions as well the total kaon cross sec- 
tions in pN and ttN collisions are of relevance. The pion 
production cross section from NN interactions is based on 
the parametrization of the experimental data by Ver West 
and Arndt |]50| and implemented in the transport model 
as described in Ref. The cross sections for the chan- 
nels ttN — !■ KY, where Y stands for a hyperon A, E, are 
taken from the analysis of Huang et al. |Q and essentially 

^ The actual final states are chosen according to the 2, 3, or 
4-body phase space. 



The question arises, furthermore, about the in-medium 
production cross sections - essentially at density po - if po- 
tentials or self energies are involved. Here we employ the 
assumption that the production matrix element squared 
|Af p does not change in the medium and the change of 
cross section can be described by a change of the avail- 
able phase space. This assumption has been employed for 
years in studies of hadron mass modifications in nucleus- 
nucleus collisions p^ . Since the experimental production 
cross sections essentially are proportional to the available 
phase space, it is sufficient to shift the threshold in a pN 
collision to 

^,^K{pN) + n'i{pA) + n'^{PK) (11) 
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using (j^), where the momenta pn ,PAtPk denote the rel- 
ative momenta with respect to the nuclear matter rest 
frame. 



2.4 Folding model and its limitations 

We briefly recall the assumptions of the folding model, 
that is used for the evaluation of momentum-dependent 
differential production probabilities from secondary pion- 
nucleon collisions in the CBUU approach and is described 
in more detail for proton-nucleus reactions in [pT||55| ,|56|. 
The underlying picture is that the proton impinging on 
a nucleus at a bombarding energy Tiab > 400 MeV is 
producing a meson x with momentum only in the first 
collision due to the available energy in the reaction. The 
Lorentz-invariant differential cross section to produce a 
meson in a primary proton-nucleon (pN) collision is given 
by H 



E. 



(fa. 



NN 



prira. 



(Tpduj E: 



(12) 

where the Pauli-blocking factor for the final nucleon states 
is neglected since kinematically the nucleons end up in an 
unoccupied regime in momentum space at the bombard- 
ing energies of interest. In eq. ( p^ S{p,u!) stands for the 
target spectral function (normalized to 1) which can be 
taken from experiment, e.g. Ref. fs^ , or parametrized ac- 
cordingly. In eq. ( p^ ) the primed indices denote meson mo- 
menta in the individual nucleon-nucleon cms frame which 
have to be Lorentz-transformed to the detection frame. 
The quantity ^/s is the invariant energy of the individual 
NN system, while the elementary differential cross section 
d^^a^ {^/s)/d^^kx in ( p^ is parametrized according to phase 
space (cf. |l^) and assumed to be isotropic in the NN cms 
frame. The momentum differential production probability 
per NN collision is obtained from (|l2|) when dividing by 
the total NN cross section. 

In order to obtain the inclusive differential kaon cross 
section in a p -I- ^ reaction within the folding model one 
has to multiply the differential cross section ( p^ ) by the 
number of first-chance collisions Ni{A). This number is 
approximately given by the area of the target divided by 
the pN cross section, i.e. Ni{A) « TTR^^j.g^f/apN- To be 
more accurate one can use Glauber theory which leads to 
Ni 7.3 for p + 12(7 ||,||. The contributions to the K+ 
yield from secondary or further sequential NN collisions 
is discarded in the folding model at subthreshold energies 
due to the energy straggling of the impinging proton and 
due to the fact that already the first-chance collisions only 
give a minor contribution to the yield observed exper- 
imentally. Note, however, that within the CBUU approach 
employed here the number of collisions is calculated dy- 
namically and thus no assumption or separate model for 
Ni{A) has to be invoked. 

Apart from the primary reaction channels described 
above, the first NN collision may also lead to the excita- 
tion of a A- resonance or even higher baryon resonances 
(e.g. N(1440), N(1535) ..) which decay to nucleons and 



essentially pions due to their short lifetimes of w 1 fm/c 
or collide with another nucleon before decaying. The AN 
or resonance-nucleon reactions are treated in the CBUU 
approach dynamically as described before. Energetic sec- 
ondary pion-nucleon collisions, however, suffer from very 
low statistics in transport simulations of p + A reactions. 
On the other hand, their contribution to the kaon yield is 
expected to be high at subthreshold energies jl^ . We thus 
implement the secondary pion-nucleon production chan- 
nels for kaons following concepts of the folding model. 

The differential cross section to produce a A'+ meson 
in a collision of an on-shell pion with a nucleon from the 
target at invariant energy then is given by 



Ek- 



N 
K 
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similar to ( 12 ) , where now the differential cross sections for 
the reactions ttN K~^Y enter. Here, the index c denotes 
all individual channels while the hyperon Y stands again 
for a yl or Z" baryon. 

In order to evaluate the yield from secondary t:N 
K^Y collisions in the folding model one folds the primary 
pion distribution - given by the primary differential pion 
cross section that is divided by the total pN cross section 
CTtot - with the nucleon spectral function S{p, uj) and the 
invariant production cross section, i.e. 



Ek 
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In Eq. (y^ the single prime indices denote the system of 
the intermediate pion and a target nucleon. Moreover, 
E'^d^apN^-n:x{"/s') / d^k'^ stands for the 7r-meson differen- 
tial cross section, which is calculated in line with (p2|), 
while atot denotes the total proton-nucleon cross section. 
The factor g^r {A) in (|^ accounts for the probability that 
the pion interacts again with a target nucleon (cf. [^5| ). 
Note, that by Eq. ( p^ the intermediate pion is assumed 
to be on-shell which might not hold for deep subthresh- 
old kaon production. As in case of the primary channel 
the expression ( |l^ ) has to be multiplied by Ni{A) in the 
folding model. Furthermore, some estimate for the sec- 
ondary rescattering probability {A) has to be employed 
in the folding model as e.g. described in Ref. j5^. As men- 
tioned before, when using the differential probabilities in 
a perturbative transport approach, the pion rescattering 
probabilities are calculated dynamically without approxi- 
mation on gT^{A). 

The one- and two-step folding model has been used 
often in the analysis of kaon or 77-meson production 
|l^,^|5^,^,^,|6^ initially employing parametrized mo- 
mentum distributions (cf. Refs. jl7[|l^) instead of spectral 
functions. They allow for an estimate of differential cross 
sections in case of weakly interacting probes and may serve 
as reference calculations for more involved simulations em- 
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ploying all initial and final state interactions during the 
reaction as the CBUU approach employed here. 



p + Pb ^ K*+X @1.5GeV 




500 1000 

Pg.+ (MeV/c) 

Fig. 3. Comparison of the differential A"*" spectra for p-\-Pb at 
1.5 GeV at different angles in the laboratory from 15° — 80° . 
The solid lines are obtained from CBUU calculations including 
kaon rescattering whereas the dashed lines are without rescat- 
tering. 



2.5 rescattering 

As indicated above, the folding model is useful for the eval- 
uation of total yields, however, becomes questionable 
in case of differential spectra especially for heavy targets 
like Au or Pb since kaon elastic rescattering cannot be de- 
scribed in a straightforward manner. In order to show the 
influence of rescatteringon the kaon spectra at differ- 
ent angles we show in Fig. y a comparison of our CBUU 
calculations with (solid lines) and without (dashed lines) 
kaon rescattering for p + Pb at 1.5 GeV. As can be seen 
from Fig. ^ the yield in forward direction {0 < 15") 
is reduced by up to a factor of 2, while for large angles in 
the laboratory (80°) the kaon spectra become enhanced 
and also shifted to lower momenta in the laboratory. Thus 



rescattering has to be included as a necessary ingredient 
for the calculations if comparisons to differential spectra 
are made on an absolute scale and especially, if one at- 
tempts to extract kaon potentials from the spectral shape 
(see below). We note that kaon rescattering thus will al- 
ways be included in the calculations to be shown in the 
following. Furthermore, we mention that the results from 
the folding model ( [l^ ) , (|lj) agree with the resulting spec- 
tra from the CBUU calculation for p + Pb and p-|~C at 1.5 
GeV within 30% when neglecting final state interactions 
as well as nuclear and Coulomb potentials. 



3 Comparison to experimental data 

In this Section we compare our calculations to the exper- 
imental spectra available from 1.0 GeV — 2.5 GeV 
bombarding energy on different targets. In order to have 
an identical assignment of lines in this section the dot- 
ted lines in Figs. 0-9 correspond to CBUU calculations 
without baryon and kaon potentials, the dashed lines show 
the results with baryon potentials included while the solid 
lines reflect calculations including both, nuclcon and kaon 
potentials as specifled in Fig. ||. In all calculations, how- 
ever, the Coulomb potential will be included by default. 

We start in Fig. ^ with a comparison to the data of the 
KaoS/SPES3 Collaboration for the differential spec- 
tra iorp+Pb andp-l-C at 1.5 GeV and 9iab = 40± 5° from 
SATURNE The experimental spectra for the Pb tar- 
get are seen to be described roughly within the error bars 
for all calculations, i.e. with and without potentials, such 
that one is not very sensitive to in-medium potentials at 
1.5 GeV in the momentum range above 350 MeV/c. For 
lower kaon momenta the repulsive potential leads to 
a sizeable decrease (or shift) in the spectra which can be 
attributed to the additional acceleration of the kaons by 
the nuclear potential when propagating from the nu- 
clear interior to the vacuum. In case of the ^^C target the 
effects from the momentum-dependent nucleon potentials 
as well as from the potential are similar to the Pb 
target. It both cases the calculations without potentials 
(dotted lines) slightly overestimate the data, but it is not 
possible to draw already some conclusion on the actual 
size of the potential since a single comparison might 
suffer from systematic errors. 

In Figs. H and we compare the CBUU calculations for 
the differential K" spectra for p + Pb and p + NaF sA 2.1 
GeV with the experimental data from the LBL Berkeley 
for different laboratory angles from 15 - 80°. At this 
bombarding energy the nucleon and A final momenta on 
average are above 0.6 GeV/c such that their potentials at 
finite density (cf. Fig. are repulsive. As a consequence 
the calculated kaon yield decreases when including the 
baryon potentials in the final state. Taking into account 
additionally the repulsive potential decreases essen- 
tially the spectrum for momenta below 250 MeV/c, but 
leaves the spectrum practically unmodified for higher mo- 
menta since the relative change of the energy by the 
kaon potential is only small. 
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Fig. 4. Comparison of the CBUU calculations for the differen- 
tial spectra for p+Pb and p+C at 1.5 GeV and 9iab = 40± 
5° with the experimental data from Ref. . The dotted lines 
are obtained from CBUU calculations without baryon and kaon 
potentials, the dashed lines show the results with baryon po- 
tentials included while the solid lines correspond to calculations 
including both, nucleon and kaon potentials. Note, that the ef- 
fect of the repulsive kaon potential is a reduction of the total 
K'^ yield as well as a shift of the spectra at lower momenta. 



In all approximations the experimental spectra are un- 
derestimated by factors ~ 2 — 3 at 15° and 35°, which at 
first sight might be attributed to an improper energy de- 
pendence of the calculations. However, a comparison of 
the CBUU calculations for p + Au with the preliminary 
data from Ref. |1 (taken at GSI) and p C at 2.5 GeV 
and Oiab — 40± 5'^^^with the experimental data from SAT- 
URNE in Fig. shows that these spectra are overes- 
timated by up to a factor of 2-3 at higher kaon momenta. 
Note, that the corresponding data from Ref. so far 
have to be considered as preliminary. These findings sug- 
gest that the spectra from Ref. |Q are systematically too 
high by about factors of 2 — 3 or the data from the KaoS 
Collaboration too low (by about the same factor). 

To further test the (over-Vunderprediction of the trans- 
port model we show in Fig. H a comparison of the CBUU 
calculations for the differential spectra for p+C at 1.2 
GeV with the experimental data from Ref. jS^l at Oiab = 
90'^ (open circles) taken at CELSIUS and Oiab = 40 ± 5° 
(full squares) from the KaoS/SPES3 Collaboration Q 
taken at SATURNE. In this particular case the spectra 
from Ref. Q are slightly overestimated by the calcula- 
tions when neglecting the kaon potential while the spectra 
from Ref. |Q at 90° are underestimated by about a factor 



p + Pb ^ K*+ X @2.1 GeV 




Fig. 5. Comparison of the CBUU calculations for the differ- 
ential spectra for p -I- Ph at 2.1 GeV with the experimen- 
tal data from Ref. |34] at different angles in the laboratory. 
The dotted lines are obtained from CBUU calculations without 
baryon and kaon potentials, the dashed lines show the results 
with baryon potentials included while the solid lines correspond 
to calculations including both, nucleon and kaon potentials. 



of 5-6 when neglecting the repulsive K'^ potential and by 
about an order of magnitude for the repulsive kaon po- 
tential included, which leads again to a sizeable decrease 
of the spectra at low momentum. It is not clear to the 
authors where such discrepancies might come from. 

The CBUU calculations demonstrate that the kaon 
spectra at 90° and 40*^ are slightly enhanced (dashed lines) 
when taking the nucleon potential effects into account. 
Contrary to the kinematical situation at Ti^b— 2.1 (or 1.5) 
GeV the nucleon and A final momenta here on average are 
below 0.6 GeV/c, where the potentials are attractive (cf. 
Fig. 0), such that the K'^ production becomes enhanced 
(dashed lines) relative to the free case (dotted lines). When 
including additionally the overall repulsive kaon potential 
the spectrum drops again (solid lines). 

We now turn to the kinematical conditions of the ANKE 
experiments at COSY-Jiilich that have taken 
spectra in forward direction for 9iab < 12°. The calcu- 
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Fig. 6. Comparison of the CBUU calculations for the differ- 
ential if"*" spectra for p + NaF at 2.1 GeV with the experi- 
mental data from Ref. |^| at different angles in the labora- 
tory. The dotted lines are obtained from CBUU calculations 
without baryon and kaon potentials, the dashed lines show the 
results with baryon potentials included while the solid lines 
correspond to calculations including both, nucleon and kaon 
potentials. 



lated differential spectra for p -|-^^ C at 1.0 GeV for 
Slab < 12° are displayed in Fig. || in comparison to the data 
from Ref. ||3^. The dotted lines again are obtained from 
CBUU calculations without baryon and kaon potentials, 
the dashed lines show the results with baryon potentials 
included while the solid lines correspond to calculations 
with both, nucleon and kaon potentials. At this low bom- 
barding energy the net attractive baryon potentials in the 
final state enhance the if + yield by about a factor of 2 
whereas the additional repulsive potential leads again 
to a decrease by a factor ^ 3. The data from Ref. |3|] are 
rather well described by the calculations that include the 
baryon and K'^ potentials (solid line), whereas the other 
limits clearly fail. This might be considered as a first indi- 
cation for the observation of a repulsive K'^ potential in 
p + A reactions, however, a full systematics in target mass 
A and laboratory energy Tiab will be needed to pin down 
this effect unambiguously. 



T„ = 2.5 GeV 
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Fig. 7. Comparison of the CBUU calculations for the differ- 
ential spectra for p + Au (upper part) and p + C (lower 
part) at 2.5 GeV and 9iab ~ 40± 5° with the experimental 
data from Refs. [59|^ . The dotted lines are obtained from 
CBUU calculations without baryon and kaon potentials, the 
dashed lines show the results with baryon potentials included 
while the solid lines correspond to calculations including both, 
nucleon and kaon potentials. 



The shaded area in Fig. 9 indicates the contributions 
from the two-step mechanisms AN K^YN and ttN — > 
K^Y, respectively, for the case of nucleon and kaon po- 
tentials included (solid line). Thus the role of secondary 
(Z\ and pion induced) reaction channels is clearly visible 
from Fig. 9 by comparing the shaded area to the solid 
line, that correspond to the total spectra. At Tiab = 1.0 
GeV the dominant fraction of the yield is due to the 
secondary channels in line with the earlier calculations in 
Refs. [T7[[l^ . Consequently, one does not probe high mo- 
mentum components of the nuclear wave function by 
spectra in a direct way. 

Without explicit representation we mention that at 
Tiab — 2.3 GeV the secondary channels in case of a Au 
target amount to about 30% and in case of a C target 
to less than 20%. This relative change with target mass 
number is attributed to the fact that for the small ^^C 
target only a fraction of the high energy pions rescatters 
in the target and produces K^Y pairs. Moreover, the role 
of the secondary channels decreases with increasing kaon 
momentum such that the high momentum tail of the 
spectra is dominated by the first chance pN production 
channel as found out before by Paryev 
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Fig. 8. Comparison of the CBUU calculations for the differ- 
ential K'^ spectra for p + C 1.2 GeV with the experimental 
data from Ref. |3q] at diab ~ 90" (open circles; lower lines) 
and diab = 40 ± S'^'^jsHj (full squares; upper lines). The dotted 
lines are obtained from CBUU calculations without baryon and 
kaon potentials, the dashed lines show the results with baryon 
potentials included while the solid lines correspond to calcula- 
tions including both, nucleon and kaon potentials. 



Furthermore, it is worth to point out that the contri- 
bution of the primary channel pN K'^YN is enhanced 
by up to a factor of 3 within the angular range of ^ < 12" 
as compared to the angle integrated yield at the energy 
Tiab = 1-0 GeV. Thus the experimental mass dependence, 
when expressed in terms of a scaling ~ A" , does not allow 
to disentangle the relative contribution of the different re- 
action channels in a satisfying manner for narrow cuts in 
the angular distribution. 

In order to provide some guideline for extrapolations 
between experiments measuring K'^ spectra at different 
angles in the laboratory we show in Fig. ^ the angular 
distribution of the kaons for momenta 0.2 GeV/c < pK < 
0.5 GeV/c as calculated within the CBUU approach for 
both, baryon and kaon potentials for p +^'^ C at Tiab = 
1.0, 1.2, 1.5, 1.8, 2.0 and 2.3 GeV. These angular distribu- 
tions are rather flat within the angular acceptance of the 
ANKE spectrometer of ^ 12*^, however, drop substantially 
for angles larger than 40^^. Thus our calculations (cf. also 
Fig. 8) do not support the idea of a 'thermal' production 
mechanism for kaons in case of p + A reactions. 




250 
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750 



Pg^+ (MeV/c) 

Fig. 9. The calculated differential spectra for p -f C at 
1.0 GeV for Oiab < 12° within the acceptance of the ANKE 
spectrometer in comparison to the data from [Q. The dotted 
line is obtained from CBUU calculations without baryon and 
kaon potentials, the dashed line shows the results with baryon 
potentials included while the solid line corresponds to calcula- 
tions including both, nucleon and kaon potentials. The shaded 
area indicates the contributions from the two-step mechanisms 
AN — > K'^YN and nN — > K'^Y, respectively, for the case of 
nucleon and kaon potentials, such that the difference to the 
solid line corresponds to the primary pN production channel. 



4 Summary 

In this work we have studied the production of mesons 
in proton-nucleus collisions from 1.0 to 2.5 GeV with re- 
spect to one-step nucleon-nucleon and two-step Z\-nucleon 
or pion-nucleon production channels on the basis of a 
coupled-channel transport approach (CBUU). We have 
included the kaon final state interactions, which are im- 
portant for heavy targets like Au or Pb, and explored 
the effects of momentum-dependent potentials for the nu- 
cleon, hyperon and kaon in the nucleus. A comparison of 
the transport calculations to the experimental K'^ spectra 
taken at LBL Berkeley, SATURNE, CELSIUS, GSI and 
COSY-Jiilich has shown that the different data sets are 
not compatible with each other. Thus no clear signal on 
in-medium potentials could be extracted from our analysis 
in comparison to experimental data so far. 

However, the detailed calculations demonstrate that 
precise and complete spectra show a substantial sensitiv- 
ity to the potentials and their momentum dependence. At 
low bombarding energies of '--^ 1.0 GeV the net attrac- 
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Fig. 10. The calculated angular distribution of the spectra 
for p-hC at 1.0, 1.2, 1.5, 1.8, 2.0 and 2.3 GeV for 0.2 <Pk < 0.5 
GeV/c. The solid lines are obtained from CBUU calculations 
including both, nucleon and kaon potentials. 



live potentials for the nucleon and the /l-hyperon in the 
final state lead to a relative enhancement of the K'^ spec- 
tra while at higher bombarding energies (~ 2 GeV) the 
baryon potentials are repulsive and thus suppress pro- 
duction relative to the free case. This phenomenon should 
be seen in the excitation function of the K'^ cross sec- 
tion when varying Tiab from 1.0 - 2.5 GeV. Furthermore, 
the shape of the spectrum for low momenta in the 
laboratory is very sensitive to both. Coulomb and nuclear 
kaon potentials, since the kaons are accelerated by both 
forces when leaving the nuclear environment and prop- 
agating to the continuum. The relative strength of this 
momentum shift in the forward K'^ spectra is propor- 
tional to the square root of the sum of both potentials, 
i.e. Ap « \/2Mk{UcouI + Uk) ■ Thus the K+ spectral 
shape at low momenta (or kinetic energies Tk) allows to 
determine the strength of the potential from exper- 
imental data in an almost model independent way espe- 
cially when comparing kaon spectra from light and heavy 
targets at the same bombarding energy |6|] as a function 
of Tk ■ Since most of the K'^ spectra measured so far have 
been taken at higher momenta in the laboratory (except 
for Ref. [Q) this finding opens up interesting perspec- 
tives for the ANKE Collaboration at COSY-Jiilich, which 
has performed a systematic study of production in 
p + A reactions down to momenta of 150 MeV/c in the 
laboratory or Tk ~ 23 MeV, respectively. 
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